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Using present global warming and paleoclimatic records from climatically sensitive regions as 
a frame of reference, we infer that global temperature changes did not exceed +0.5°C during 
the current interglacial or +2.0°C during the last glacial period. In order to completely explain 
such fluctuations by solar irradiance changes, all solar and terrestrial factors must be optimized. 
Since this is unlikely, we conclude that solar forcing of pre-anthropogenic climate change is a 
significant and perhaps dominant factor but other processes must also be significant. Solar 
irradiance changes alone cannot result in the temperature decrease required to change from 
interglacial to glacial as suggested by Opik (1965). 


1. Introduction 
1.1. Glacial and interglacial climate 


During the Pleistocene (1.8 Ma to present), the Earth’s climate has varied between two 
states, glacial and interglacial. This can be seen in the oxygén isotope stratigraphy 
in Figure 1. The meteorological cycle has been likened to a multiple-stage still with 
6 to 7 theoretical distillation plates (Epstein 1958). The amount of fractionation in 
the hemispheric distillation process is a function of altitude, latitude and temperature 
gradient from equatorial to polar regions. The variations of the 180/160 ratio during 
the hydrologic cycle are measured in %o relative to mean seawater. With reference to 
the Byrd core, Antarctica, it can be seen that 6180 fluctuates around —40%o from about 
110 ka to 18 ka, gradually increasing after 18 ka to fluctuate around a new steady-state, 
—34%o, from about 10 ka to present. The increase in 6180 from —40%o to ~34%o is the 
result of global warming and coincident with rising sea-level as the continental ice masses 
melted (Fairbanks 1989). The Camp Century, Greenland core shows 6'8O variation 
around the same two steady states but during the transition form 18 ka to 10 ka, the 
warming is temporarily interrupted by a return to glacial conditions during the Younger 
Dryas from about 12,900 BP to 11,600 BP. Thus, the hydrologic cycle is responding 
to both regional and global climate changes. The North Atlantic is more sensitive to 
climate change during the glacial period as can be seen by the larger fluctuations of 
6'8O in Greenland ice cores when compared to Antarctic ice cores but the two states are 
global in nature. In general, global climate during the Pleistocene can be modeled as a 
set of systems fluctuating around the glacial and interglacial steady states separated by 
transitional episodes. 6!*O fluctuations around steady state are three to four times less 
during the interglacial compared to the glacial epoch. 


1.2. Regional vs. global climate change 


According to estimates from the CLIMAP project (CLIMAP Project Members, 1981), 
the Laurentides ice mass reached a thickness of 3800 m and the northwestern European 
ice mass reached a thickness of 2500 m, whereas the Greenland ice cap reached an 
intermediate thickness of about 3200 m. Using an empirical equation from Dansgaard 
(1964), we can estimate the glacial age temperature at Camp Century to have been 
= —38°C. The Northern Hemisphere continental areas covered by ice must have been 
at similar temperatures, warming to their present annual temperatures. As a robust 
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FIGURE 1. Ice core 6'50 (Camp Century, Greenland, after Dansgaard et al. 1982 (la); Byrd 
Station, Antarctica, after Johnsen et al. 1972, (1b)). Note that both isotopic records show . 
distinct modes between glacial and interglacial states. For the glacial state (20 ka-13 ka), 
6150 varies by an amplitude of 2%o on a century time-scale about an average value of + —40%o. 
During the interglacial (9 ka-present), 6150 varies by an amplitude of <0.5%o around an average 
of © —29%o for Camp Century and % —33%o for Byrd Station. The decreased 150 variability 
suggests a decrease in sensitivity to climate forcings. Camp Century shows greater glacial to 
interglacial contrast due the influence of surface ice and snow expanse during the glacial period. 
A sharp increase in the Camp Century record contrasts with the more gradual increase in 
the Byrd Station record. This may reflect the ameliorating influence of the Southern Ocean. 
The pronounced Younger Dryas negative deflection (10.5 ka) present in the Camp Century 
may correlate with a much smaller feature (+10 ka) in the Byrd Station record (calibrated 
ages are about 1700 years older as in the text). The Younger Dryas event may have affected 
North Atlantic circulation greatest and thus adjacent continental indicators have more profound 
signatures. Thus, 6180 records contain both global and region climate signals. 
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Region AT (°C) Area of Globe (%) Contribution to Global 
Mean weighted by area (°C) 

Presently Glaciated —10 3 —0.3 

Not glaciated at present —40 5.7 —2.3 

Other continental areas —6 20.5 —1.2 

The world oceans —2.3 70.8 —1.6 

Hansen et al., 1984 2 = —5.4 

GCM Model: x = —5.3 


TABLE 1. Temperature changes (AT) from present to 18 ka. 


but not rigorous estimate of the temperature increase for those areas, we suggest about 
40°C. For the areas still covered by polar ice, we estimate a change of about 10°C. From 
a paleoclimatology summary by Crowley & North (1991), we take the average change at 
mid-latitudes to be between 5-8°C and 4-5°C for tropical continent and areas under the 

` influence of maritime air masses. Let us assume an average of ~6°C for continental areas 
not covered by ice sheets. CLIMAP Project (McIntyre et al. 1976) estimated an average 
change for the global oceans of only 2.3°C but with great variation between different 
ocean basins. For example, in the Northern Hemisphere, they estimated temperature 
changes varied from 0.8°C for the Indian Ocean to 3.8°C for the North Atlantic. We 
summarize and weight these changes by areas in Table 1, arriving at a global change 
of —5.4°C from the glacial to the present state. This compares favorably with more 
rigorous GCM (global circulation model) modeling by Hansen et al. (1984) and serves 
to illustrate several pertinent points. First, the global oceans serve as a great source 
and sink of heat moderating climate. Second, 8.7% of the area of the Earth accounts 
for ~50% of the global warming during the transition from glacial to interglacial states. 
Third, there were large regional differences in temperature changes around the global 
mean. 

As previously noted, oxygen isotope stratigraphy demonstrates that the hemispheric 
temperature fluctuations during the past 10,000 years (the Holocene) are relatively small. 
The linear empirical relationship between 6180 and temperature (Bradley 1985) suggests 
that these fluctuations do not exceed £0.6°C (or 1.2°C peak-to-peak). If this is correct, 
global warming during the past 100 years of 0.5°C based on instrumental records (Folland 
et al. 1990) should be a very significant event relative to previous events during the 
Holocene. Of course, regional temperature variations of sub-hemispheric scale have been 
both lesser and greater. For example, there was an insignificant increase in temperature 
during the last century for the contiguous United States comprising roughly 1.6% of 
global area (Hanson et al. 1989). However, from data in Trends, 1991 (Boden et al. 
1991), we calculate an increase of (1.2+0.2)°C(s) from solar cycle 14 (1901-1913) to 
solar cycle 21 (1976-1987) for the southern deserts which is one of 23 regions of the 
U.S.A. reported in Trends, 1991. These sorts of regional variations of climate change 
together with conceptual and dating problems make it difficult to evaluate past changes 
in Holocene global temperatures. The instrumental record coupled with GCM studies 
can be helpful in resolving this problem. Specifically, in order to evaluate the feasibility 
of solar forcing of global climate, we need to know the magnitude of global Holocene 
temperature changes to be forced. As noted above, oxygen isotope studies of polar ice 
suggest global changes equal to or less than +0.6°C during the Holocene. This suggestion 
needs further justification. 

GCM studies suggest that global temperature change will be greatest at high altitudes 
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and high continental latitudes (Wang et al. 1991). High altitude temperatures are 
governed by evaporation of surface ocean water with transfer of 540 calories per gram 
of latent heat to higher latitudes. Albedo changes, as the result of changes in the ice 
and snow cover, produce high-latitude continental sensitivity to global climate forcing. 
Evidence for climatic events like the Medieval Warm Epoch and Little Ice Age come 
from North Atlantic border-lands and Europe above 45°N. It has been possible recently to 
demonstrate by study of montane glaciers that these climatic events are global rather than 
merely regional in scope (Rothlisberger 1986). Montane glaciers also provide perspective 
concerning the relative intensity of 20th Century global warming. 

Hasenrath & Kruss (1992) have compared glacial maps for Mount Kenya (elevation 
4.6 to 5 km). They estimate that the Mount Kenya Glacier retreated by 40% in area 
from 1963 to 1987. The retreat of the terminus of the Speke Glacier in the Ruwenzori 
Range of Uganda has accelerated from an average of +2 m/yr between 1958 and 1977 to 
more than 11 m/yr between 1977 and 1990 (Kaser & Noggler 1991). The areal extent of 
worldwide montane glaciation has been reduced by about 50% during the last century 
(Thompson 1991). Thompson et al. (1989) have studied the Dunde Ice Cap on the 
Qinghai-Tibetan Plateau at a summit elevation of 5235 m. According to these authors, 
in that region “the last 60 years were apparently one of the warmest periods in the entire 
record, equaling levels of the Holocene maximum between 6000 and 8000 years ago” (p. 
474), 

Evidence for high-latitude warming comes from dendroclimatologic studies and mea- 
surements of the temperature gradients in bore holes. Innes (1991) has documented high- 
altitude and high-latitude increases in 20th Century tree growth from about 30 studies of 
trees from Northern Europe and North America. In contrast to Thompson’s conclusion, 
Innes states, “It appears that neither the northern extent nor the present growth rates of 
boreal forests have reached the extent of rates characteristic of the mid-Holocene, when 
the climate was demonstrably warmer.” However, from their dendroclimatologic study 
of high-latitude boreal sites in Canada and Alaska, Jacoby & Arrigo (1989) conclude 
that “the global warming trend over the past century of increasing CO% has exceeded 
the recent (past three centuries) level of natural variability (pre-anthropogenic) of the 
climate system”. In a study of old-growth Huon pines from a sub-alpine site in western 
Tasmania, Cook et al. (1992) concluded that, “Although the most recent 25-year period 
is warmer than any comparable period over the past 1090 years, it is not yet statisti- 
cally unprecedented” (p. 215). Graybill & Shiyatov (1992) have reconstructed June-July 
temperatures for ad 961-1969 using tree-ring width variation of Siberian Larch from the . 
Polar Urals (65-67°N). They concluded that “the relatively sustained rise in tempera- 
ture from the 1880s to the 1960s has the greatest magnitude in the record for such an 
extended period (1.6°C difference between the 20-year intervals centered on 1870 and 
1960)” (p. 411). Note that for this climate sensitive area, this recent regional warm- 
ing event is threefold greater than 20th Century global warming of ~0.5°C. Bore hole 
temperature gradients from Central and Eastern Canada (Beltrami et al. 1992; Shen & 
Beck 1992; Wang et al. 1992) as well as from the Alaskan Arctic (Lachenbruch & Vaughn 
1986; Kakuta 1992) confirm the temperature rise from the late 19th Century through 
the late 20th Century. The temperature rise in Canada was typically 1-3°C and greater 
at Prudhoe Bay, Alaska where Kakuta (1992) estimates warming by 4°C from ad 1853 
to ad 1973. 

Summarizing, the modern instrumental record has made it possible to estimate 20th 
Century global warming of about 0.5°C (Folland et al. 1990). There is insufficient data 
to establish directly global climate change prior to the 19th Century. The existing data is 
biased toward relatively high-latitude, climatically-sensitive regions of the North Atlantic. 
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However, our comparison of the recent warming event with the paleoclimate record within 
climatically-sensitive regions suggests that this recent event exceeds in intensity natural 
global climate variability during the last millennium and is of comparable magnitude to 
the Holocene maximum between 6000 and 8000 years ago. Oxygen isotope stratigraphy 
also demonstrates the relative global stability of pre-anthropogenic climate. Solar forcing 
needs only account for global temperature variations of no more than +0.5°C to account 
for all the natural climate variability during the last millennium. The Holocene maximum 
may exceed a global temperature increase of +0.5°C. However, our previous estimate of 
+0.6°C variability from the polar oxygen isotope record is probably a conservative upper 
limit for global temperature change during the Holocene. 


2. Solar forcing of Holocene climate 
2.1. The Sun as a low-frequency harmonic oscillator 


The production of cosmogenic 14C is modulated by solar activity via the magnetic field 

` associated with the solar wind. The relationship is inverse and, so, the various histori- 
cally observed solar minima such as the Maunder Minimum occur as *C maxima and 
the Medieval Solar Maximum from ad 1100-1250 occurs as a '*C minimum in the spec- 
trum of !*C variations. The reader may refer to two papers for recent reviews of the 
subject (Damon & Sonett 1991; Stuiver et al. 1991). Damon & Sonett (1991) concluded 
that medium- and short-term variations with periods of ~208 (Suess), 88 (Gleissberg), 
and 11 (Schwabe) years are in large part produced by solar wind modulation. A longer 
period of +2300 years (Hallstattzeit) seemed enigmatic but appears to 100% amplitude 
modulate the Suess cycle. The ~2300-year period seems to be associated with recurrent 
Little Ice Ages (Schmidt & Gruhle 1988). Stuiver & Braziunas (1991), using the maxi- 
mum entropy spectral method (MEM), chose a 420-year period as a fundamental with 
harmonics related to changes in the solar convective zone and noted that it splits into 
three periods at higher AR order. We note that the apparent splitting is an artifact of 
the MEM method. The peaks are always present but vary greatly with the autoregres- 
sive (AR) order. Too low an AR order results in oversmoothing. We deduced from both 
DFT and MEM power spectral analyses that the four periods (526, 424, 356 and 300 
years) are the 4, 5, 6, 7 harmonics of the poorly resolved Hallstattzeit period that we 
estimated from its harmonics to be 2120 years (Damon & Jirikowic 1992a,b). Further 
analysis yielded 212 years for the Suess period with strong (>95% confidence) harmonics 
at 106 years (2nd harmonic) and 53 years (4th harmonic). 

The harmonics of the Suess period combined with the Gleissberg period appear to 
strongly modulate the 11.1-year (Schwabe) cycle. Furthermore, a reasonable fit to the 
annual Wolf Sunspot Indices can be obtained using a squared amplitude modulation 
model (after Sonett 1982) by allowing these periods (212, 106, 53, 88.3 years) to amplitude 
modulate the 11.1-year carrier (Damon & Jirikowic 1992a). 

Subsequent to the work summarized above, one of the authors (Jirikowic) added Thom- 
son’s (1990) multi-taper power spectral analysis technique to the techniques made avail- 
able to us by Prof. Charles Sonett here at the University of Arizona (periodogram, 
Bretthorst-Bayesian, Maximum Entropy [MEM], Discrete Fourier Transform [DFT)]). 
The analysis of the latest version of the high-precision 14C data set (Stuiver & Reimer 
1993) is shown in Figure 2. The Hallstattzeit fundamental and first nine overtones (2- 
10fy) clearly resolve at >95% confidence (harmonic F-test). The Suess period (fg and 
10fy) has unusually high power spectral density and appears to act as a fundamental 
as well as a harmonic of the Hallstattzeit. The accounts for the high spectral power 
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FIGURE 2. Power spectral density estimates using Thomson’s Multi-taper Method (Thom- 
son 1990) and three autoregressive modeling orders: 120 (dotted), 160 (solid) and 200 (dashed) 
from the 1993 Radiocarbon Calibration tree-ring 5'*C data (Stuiver & Reimer 1993). Annotated 
peaks show assignment to hypothetical harmonics to the fundamental Hallstattzeit (subscript 
H), Suess (subscript S) and Gleissberg (subscript G) variations. All these peaks have greater 
than 95% F-test confidence. Although the Hallstattzeit is poorly resolved, the annotated har- 
monics combine to suggest a fundamental period of 2108 + 15 years. Earlier spectral analyses 
(Damon & Jirikowic 1992b) showed a strong fourth harmonic of the Suess that is not apparent 
here. 


spectral density of its second harmonic (2fs, >95% confidence). However, its fourth har- 
monic appears as a relatively weak doublet, unlike the MEM and DFT techniques (>95% | 
confidence). The Gleissberg period is prominent at >95% confidence (88.1 years). 

Neglecting the poorly resolved second Hallstattzeit harmonic but using harmonics 3-10 
plus 14, 15, 17, 18 (all with >95% confidence) to calculate the Hallstattzeit fundamental, 
we obtain 2108+15 (ő years compared with 2115+18 (@) years previously obtained using 
the combined MEM and DFT methods (Damon & Jirikowic 1992b). 

The enigmatic Hallstattzeit period may be divided into two episodes. For +1520 
years, it is associated with high solar activity and weak periodicity followed by a ~600- 
year episode in which solar activity is low on the average and the Suess period and 
its harmonics as well as the Gleissberg period are strong. These 600-year episodes are 
associated with Little Ice Ages. Thus, when the Medieval Solar Maximum ended at 
ad 1250, the Little Ice Age began accompanied by a series of profound solar minima 
and terminated after about 600 years toward the middle of the 19th Century when the 
Modern Solar Maximum began (Jirikowic & Damon 1994). 

Instrumental analysis has confirmed a significant change in solar irradiance of 0.1% 
during the 21st solar cycle (Willson & Hudson 1988, 1991). Paleoclimate data suggests 
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Change in A Forcing ATeq (°C) ATeq (°C) ATeq (°C) 


Solar Solar’ (tropopause Glacial? Present Glacial? Future? 
Condition Irradiance level) 2.4°C 1.0°C 0.6°C 
Present Sun 0.0 0.0 0.0 0.0 0.0 
AD 1980-1986 

No plage, spots —1.5 —0.36 —0.9 —0.4 —0.2 
or network 

No-cycling stars —2.4 —0.58 —1.4 —0.6 —0.3 
(average) 

Minimum possible —3.5 —0.84 —2.0 —0.8 —0.5 
Sun 


"From Lean et al. (1992). 
2Calculated from Hansen et al.(1984). 
3Calculated from Hansen et al.(1984) assuming variations are negligible on an ice-free Earth. 


TABLE 2. Decreasing equilibrium temperature (T.,) as a function of solar condition. 


that longer periods such as the Hallstattzeit, Suess and Gleissberg periods may also be 
associated with changes in solar irradiance (e.g., see Eddy 1977). Further evidence comes 
from solar and astrophysics. 


2.2. The Sun as a variable star 


As discussed elsewhere in this colloquium’s proceedings (see Lockwood 1994), 74 solar- 
type stars have been observed at the Mount Wilson Observatory. Of these 74, 13 have 
been studied monthly since 1966 and nightly since 1980 (Baliunas & Jastrow, 1990). Only 
sporadic measurements were made since 1978 for the other 61 stars. Nine of the 13 stars 
were observed to show a range of magnetic activity similar to the solar cycle whereas four 
that showed no variation of magnetic activity with time were assumed to be analogous 
to the Sun during the Maunder and other solar minimum. The nine magnetically active 
stars had regular activity cycles like the Sun. Their brightness varied by up to 0.4% 
compared to 0.1% for the Sun during cycle 21. Baliunas & Jastrow (1990) state, “Our 
results suggest that when a solar-type star enters or leaves a Maunder Minimum, its 
brightness may change by more that the 0.1% observed during the Sun’s 11-year cycle” 
(p. 522). Lockwood et al. (1992) presented a compilation of eight years of observations 
of 33 stars at the Lowell Observatory. They report that the stars whose mean activities 
were closest to the Sun’s varied in brightness by 2.5 to 10 times the variation in brightness 
of the Sun. They state, “This raises the question of whether the Sun’s cycle irradiance 
variation is always so small or whether larger changes occur in other cycles or over longer 
time intervals” (p. 655). 


2.3. Solar variability and the Maunder Minimum 


Based on current knowledge of the Sun and Sun-like stars, Lean et al. (1992) estimated 
the change in solar irradiance and forcing at the top of the atmosphere for various 
conditions of the Sun (see Table 2). Assuming, as did Baliunas & Jastrow (1990), that the 
Maunder Minimum was equivalent in change of irradiance to the average solar-like non- 
cycling stars yields a 2.4 change in irradiance. Also, assuming that a global equilibrium 
temperature decrease of 1°C was required to explain Little Ice Age cooling, they conclude, 
“.. that solar forcing may have contributed some, but not all, of the observed warming 
from the Little Ice Age to the present time” (p. 1594). 

Keeping in mind that the ~0.5°C rise in temperature in recent time began from a 
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Feedback Gain AT (°C) 
Unity 0 -0.6 
Water vapor—clouds 0.415 —2.2 
Land ice 0.169 —0.9 
Sea ice 0.132 —0.7 
Vegetation 0.057 —0.3 
CO2 0.113 —0.6 


gi = 0.886 SO ATeq = -5.3 


TABLE 3. Gains (gi) and temperature increments resulting from feedbacks at the Ice Age 
Maximum (18 ka). 


relatively cool global temperature in the late 19th Century, and this rise approaches or 
exceeds previous such warming events during the Holocene (see 1.2 above), we suggest 
that the estimate of 1°C decrease in global cooling exceeds what is required to explain 
the Little Ice Age cooling during the Maunder Minimum. We concluded that less than 
+0.5°C global cooling is required to explain climate variability during the last millen- 
nium. The global temperature rise from the Maunder Minimum to ad 1980-1986, based 
on the paleoclimatic record and instrumental measurements discussed, is greater than 
0.5°C but less than 1°C. 

In order to evaluate the temperature decrease resulting from a given change in solar 
irradiance (Table 2), various positive feedbacks must be considered. For our analyses, 
we refer to Hansen et al. (1984). Their paper included analysis of the various gains (g;) 
and temperature increments (AT;) including those for water vapor-clouds, land ice, sea 
ice, vegetation, and CO% (see Table 3). 

From Table 3, the equilibrium temperature change, AT.g=5.3°C and the temperature 
increment that would occur without feedbacks, ATp=0.6°C. The total feedbacks, S* f; 
is related to the total gain, J` g:, by: 


1 | 
oe I oss 88 (2.1) 


t 


and AT is related to AT, by: 


AT, 33° 
ATp = pee 60 (2.2) 


eq = 
5 fi 8.8 
The relative increment in solar energy required to yield the increment ATbp is: 
ASo ATp 4 x 0.6°K -3 
S i æ ie) 


Note that the surface temperature of Earth is warmer by ~33°K due primarily to the 
greenhouse gases. The change in solar energy required to obtain ATp is 


9.4 x 107? x 1367.54Wm~? = 12.8Wm~” 
that corresponds to a forcing at the tropopause level of 12.8 Wm~?+5.7 or 2.2 Wm~? 
This yields a climate sensitivity of 5.3/2.2 Wm-? = 2.4°C/Wm7?. 
Likewise, Hansen et al. (1984) estimate > f; at present to be 3 to 4, we take f ~3.6 


that yields a climate sensitivity of 1°C/Wm~?. If we assume in a world without ice, 
albedo effects are negligible, the temperature increase required to attain an ice-free state 
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f __—— Recot Nie 


a ice Free Earth ——__—-> Maaah 


f=positive feedback factor 


FIGURE 3. The relation between equilibrium surface temperature (Teq), positive feedbacks 
(© fi) and climate sensitivity. As a consequence of the increased positive feedbacks in colder 
climates, the same change in solar energy would have the greatest effect during glacial epochs 
and the least effect for an ice-free Earth. 


would be a minimum of 4.8°C (Barron & Washington 1984), the sum of feedbacks would 
be 2.1 and the climate sensitivity would be ~0.6°C/Wm~?. 

In order to go from the present climate to an ice-free (swamp-ocean) climate requires a 
forcing at the tropopause level of 8.6/Wm~? or an increase of 3.5% in the solar constant 
whereas a cooling from present to the climax of the last ice age at 18 ka requires only 
a forcing decrease of 2.2 Wm~? or a decrease in the solar constant of 9.4x10-3 %o. 
Neither of those changes are obtainable by changes in solar irradiance alone (Table 2). 
However, significant solar modulation of climate change within the two steady states, 
glacial-inter-glacial, can be obtained. Neglecting the increased positive feedback during 
a Little Ice Age, which we estimate to be only ~10%, a peak-to-peak change of 0.6°C 
(Table 2) explains most, if not all, of the Holocene temperature variation. The same 
change in solar irradiance, 2.4%, would produce a peak-to-peak variation of 1.4°C at 
18 ka. Figure 3 illustrates the increased variability of climate during the ice ages and 
decreased variability at present and during an ice-free climate. 

The evidence that the recent warming event approaches or perhaps exceeds past events 
sets limits on natural variability during the Holocene but we must keep in mind that 
warming since the Little Ice Age is the combined effects of increased solar activity and 
greenhouse forcing. According to Hansen & Lacis (1990), greenhouse forcing to 1989 was 
2.04 Wm7?, to which must be added our estimated solar forcing of 0.58 Wm~? (Table 2) 
or about 2.6 Wm~? total. This forcing at equilibrium, assuming the validity of the climate 
sensitivity of 1°C/Wm~?, would yield an increase of 2.6°C since the Maunder Minimum. 
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However, heat is transported to the mixed layer of the ocean and from there to the deep 
sea resulting in a long lag time. From their three-dimensional (3D) model, Hansen et 
al. (1984) estimated “the time for the global area-weighted mixed layer temperature 
to reach 63 percent of its equilibrium response in 124 years” (p. 156). The 124-year 
e-folding time is applicable for the high climate sensitivity of 1°C/Wm~?. If we divide 
the 0.6°C (Table 2) equilibrium temperature increase equally between warming after 
the Maunder Minimum and additional warming during the Modern Solar Maximum, we 
calculate that 25% of the current warming of 0.5°C is from solar forcing and 75% is the 
result of greenhouse gas forcing. The temperature rise following the Maunder Minimum 
would attain its equilibrium value of 0.3°C to which 0.125°C must be added from the 
warming due to the current solar maximum. The resultant 0.425°C would then constitute 
a little more than half of the 0.8°C warming since the Maunder Minimum. 

Our argument, thus far, has been based upon the high climate sensitivity of 1°C/Wm7?. 
As a result of uncertainties in the feedbacks, the range of possible climate sensitivities lies 
between 0.5°C/Wm-~? and 1.25°C/Wm~? (Hansen & Lacis 1990). The e-folding time 
is also a function of the feedbacks (J, f;). For the lowest value of climate sensitivity, 
0.5°C/Wm~?, the e-folding time is only about 40 years. This much shorter residence 
time slightly more than compensates for the lower greenhouse forcing climate sensitivity 
but the solar contribution falls to 13% of the current global warming. The total tem- 
perature rise since the Maunder Minimum becomes 0.65°C, of which only 33% is due to 
solar forcing. 

It is of interest to apply this analysis to pre-anthropogenic climate change. For this, 
we take predicted temperature change between the Medieval Solar Maximum (ad 1100- 
1250) and the immediately following Wolf Minimum (ad 1250-1350). We assume, as 
above, that the equilibrium warming and subsequent cooling would each be 0.3°C. The 
predicted temperature fall must be calculated from the e-folding time of 124 years. It 
would be 0.4°C. This global change is not great but would be observable in climatically 
sensitive regions using appropriate paleoclimate indicators. It would probably be below 
detection levels in other regions. If the climate sensitivity is 0.5°C/Wm~?, the tempera- 
ture decrease would only be ~0.3°C but probably still observable. At high latitudes and 
at high elevations, the temperature might change by 1°C as compared to a 2°C warming 
during the current global warming event. There remains the possibility of greater cooling 
during prolonged periods of intense solar minima such as the Little Ice Ages (ad 1250- 
1850), for example. To completely account for Holocene temperature changes, however, 
it seems that all permissible solar and terrestrial factors must be maximized. 


3. Conclusions and summary remarks 


1. Pleistocene climate alternates between two steady states, glacial and interglacial, 
separated by transitional episodes. Climatic fluctuations around steady state are greater 
during glacial epochs than interglacial epochs by about a factor of three to four as 
observed in oxygen isotope stratigraphy. 

2. Regional temperature changes are much more variable than and should not be 
confused with global temperature change. Since climatic events like the Little Ice Age, 
Medieval Warm Epoch and Altithermal were originally defined based upon paleoclimatic 
and historical events from Northern Europe and the North Atlantic regions, there has 
been a tendency to exaggerate global temperature changes during the Holocene. 

3. Present global warming can serve as a frame of reference in evaluating Holocene 
global temperature changes from paleoclimate records within climatically sensitive re- 
gions such as high-latitude and high-altitude continental areas. Such comparison strongly 
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suggests that global temperature change did not exceed +0.5°C during the last 10 ka or 
+2.0°C during the last glacial period. 

4. The Sun is a low frequency harmonic oscillator with solar activity fundamental 
frequencies of 1/2120 yr, 1/212 yr, 1/87 yr and 1/11 yr plus harmonics. The 1/2120 
yr cycle includes two phases, episodes lasting 600 years during which the century-scale 
variation includes pronounced minima such as the Maunder Minimum separated by 1500- 
year episodes of high and more chaotic activity. 

5. The Sun in one of a family of Sun-like stars most of which have cycles of brightness 
similar to the Sun’s. These Sun-like stars vary in brightness up to four times the variation 
in brightness of the Sun. Some of the Sun-like stars are in quiescent states similar to the 
Sun during the Maunder Minimum (Baliunas & Jastrow 1990; Lockwood et al. 1992). 

6. Lean et al. (1992) suggest that, during a solar minimum, the Sun decreases in 
brightness by 2.4%o relative to the present Sun similar to the average variation of the 
non-cycling stars (Table 2). 

7. Based on the 3D GCM analysis of Hansen et al. (1984), decreases in irradiance 
of 2.4% during solar minima relative to the present Sun would result in decrease of 
equilibrium global temperature by from 0.3°C to 0.6°C during the present interglacial 
and 0.9°C to 1.4°C during the last glacial maximum at 18 ka. 

8. Applying paleoclimatic restraints on global temperature and using the 3D GCM 
analysis of Hansen et al. (1984) with a climate sensitivity of 1°C/Wm~?, we estimate 
that global temperature increased by ~0.3°C following the Maunder Minimum and more 
than 0.1°C of the 0.5°C of the current global warming since the late 19th Century would 
be the result of increased solar irradiance during the Modern Solar Maximum. Thus, 
about 50% of warming since the Maunder Minimum would be the result of increase solar 
activity. If the lowest estimated climate sensitivity is assumed (0.5°C/Wm 7), the solar 
contribution would be only about 33%. 

9. A climate sensitivity of 1°C/Wm~? requires that, if no more greenhouse gases were 
released to the atmosphere, global average temperature would continue to increase by 
2.6°C due 20% to solar activity and ~80% to the effect of greenhouse gases. Such an 
increase would be unprecedented during the last 2 Ma. With inevitable further release 
of greenhouse gases, global climate would tend towards an ice-free globe held in check 
only be the large thermal inertia of the global oceans and polar ice. 

10. Assuming the validity of the astrophysics-solar physics analysis of Lean et al. 
(1992), changes in solar irradiance cannot result in a change from interglacial to glacial 
as suggested by Opik (1965). However, the century-scale and millennial changes in 
solar activity can account for much of the climate fluctuations around the glacial and 
interglacial steady states due to the existence of strong positive feedbacks. 

11. In order to completely explain Holocene global temperature changes by solar 
forcing, all terrestrial and solar factors must be optimized. Since this is unlikely, we 
conclude that solar forcing of pre-anthropogenic global climate change is a significant 
and perhaps dominant factor, but other processes must also be significant. 

12. Added note: The Vostok Antarctica core demonstrates relative climate stability 
during the present interglacial (Holocene) and also the penultimate (Eemian) interglacial 
(Dansgaard et al. 1993). The apparent instability during the Eemian suggested by the 
recent GRIP and GISP cores from Greenland (ibid.) appears to be caused by layer 
thinning and deformation rather than climate instability (e.g., Grootes et al. 1993). 
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